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Abstract —Hardware and software architecture of a novel 
spacecraft Attitude Determination and Control System (ADCS) 
based on smartphones using Variable Speed Control Moment 
Gyroscope (VSCMG) as actuator is proposed here. A space¬ 
craft ground simulator testbed for Hardware-ln-the-loop (HIL) 
attitude estimation and control with VSCMG is also described. 
The sensor breakouts with independent micro-controller units 
are used in the conventional ADCS units, which are replaced by 
a single integrated off-the-shelf smartphone. On-board sensing, 
data acquisition, data uplink/downlink, state estimation and 
real-time feedback control objectives can be performed using 
this novel spacecraft ADCS. The attitude control and attitude 
determination (estimation) schemes have appeared in prior pub¬ 
lications, but are presented in brief here. Experimental results 
from running the attitude estimation (filtering) scheme with 
the “onboard” sensors of the smartphone in the HIL simulator 
are given. These results, obtained in the Spacecraft Guidance, 
Navigation and Control Laboratory at NMSU, demonstrate the 
excellent performance of this estimation scheme with the noisy 
raw data from the smartphone sensors. 

1. INTRODUCTION 

Spacecraft Attitude Determination and Control System 
(ADCS) is considered to be the vital sub-system of all other 
spacecraft systems, as it enables the critical capability of 
attitude maneuvers such as pointing, re-orientation, stabiliza¬ 
tion etc. An ADCS is necessary for commercial satellites, 
as well as spacecraft with science and tactical missions like 
observation of solar system bodies, stars and galaxies, as well 
as features and objects on the Earth’s surface. The ADCS 
system consists of the attitude determination and estimation 
component, and the attitude control system [1]. A novel 
mechatronics approach is proposed for spacecraft ADCS that 
uses a Commercial Off-The-Shelf (COTS) smartphone with 
integrated sensors as on-board computer. This technique uti¬ 
lizes the smartphone’s inbuilt accelerometer, magnetometer 
and gyroscope as an Inertial Measurement Unit (IMU) for 
attitude determination and a Variable Speed Control Moment 
Gyroscope array for attitude control. 

The primary motivation for using an open source smart¬ 
phone is to create a cost-effective, generic platform for 
spacecraft attitude determination and control, while not sacri¬ 
ficing on performance and fidelity. The PhoneSat mission of 
NASA’s Ames Research Center demonstrated the application 
of COTS smartphones as the satellite’s onboard computer 
with its sensors being used for attitude determination and its 
camera for Earth observation [2]. More recent satellites in the 
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Phonesat series, Phonesat 2.4 and 2.5, successfully utilized 
smartphones for attitude control using reaction wheels [3]. 
University of Surrey’s Surrey Space Centre (SSC) and Surrey 
Satellite Technology (SSTL) developed STRaND-1, a 3U 
CubeSat containing a smartphone payload [4], [5]. 

Eor a typical mechatronics system, a number of sensor 
breakouts including accelerometer, rate gyros, magnetome¬ 
ter, communication breakouts and GPS units are used as 
independent sensor breakouts. Signals from these breakouts 
are processed through an on-board microcontroller, which 
estimates the states and calculates the control signals based 
on the feedback control law, which are then sent to the actu¬ 
ators. A more complex mechanical system like a spacecraft 
with internal actuators demands a microcontroller with more 
processing power, memory capacity and input/output ports, 
which can considerably increase the cost and mass budget of 
the system. Evolution of smartphones have been remarkable 
in the past few years with sophisticated processors, higher 
order system architecture, and high quality sensors integrated 
into a compact unit. These smartphones can be used as on¬ 
board CPU for performing HIL operations of more complex 
mechatronic systems like spacecraft with internal actuators. 
Some advantages of using smartphones, on-board are: 

1) compact form factor with powerful CPU, GPU etc., 

2) integrated sensors and data communication options, 

3) long lasting batteries: reduces total mass budget, 

4) cheap price and open source software development kit. 

The novel ADCS scheme developed here includes an 
attitude and angular velocity estimation scheme that is based 
on inertial directions and angular velocity of the spacecraft 
measured by sensors in the body-fixed frame of the smart¬ 
phone. The state estimation scheme presented here has the 
following important properties: (1) the attitude is represented 
globally over the configuration space of rigid body attitude 
motion without using local coordinates or quaternions; (2) 
the scheme developed does not assume any statistics (Gaus¬ 
sian or otherwise) on the measurement noise; (3) no knowl¬ 
edge of the attitude dynamics model is assumed; and (4) the 
continuous and discrete-time filtering schemes presented here 
are obtained by applying the Lagrange-d’Alembert principle 
or its discretization [6] to a Lagrangian function that depends 
on state estimate errors obtained from vector measurements 
for attitude and angular velocity measurements. 

The spacecraft’s attitude is controlled by the variable 
speed control moment gyroscope (VSCMG) actuators. A 
general dynamics model of a spacecraft with VSCMGs is 
derived in [7], [8], and is adopted here. Since the configu- 



ration space of attitude motion of a spacecraft with internal 
actuators is a nonlinear manifold, the global dynamics of 
this system is treated using the formulation of geometric 
mechanics [9], [10]. This model is obtained using variational 
mechanics, and in the framework of geometric mechanics 
on the nonlinear state space of this system. This dynamics 
model relaxes simplifying assumptions used in most of the 
existing literature on (VS)CMGs [11]-[14]. These simplify¬ 
ing assumptions, which are relaxed here, includes; 

1) Center of Mass (CoM) of gimbal is perfectly aligned 
with rotor CoM along the gimbal axis {a = 0), 

2) axisymmetric rotor, 

3) both the gimbal & the rotor-fixed coordinate frames 
are their corresponding principal axes frames, 

4) rotor & gimbal inertias are about their respective CoM, 

5) gimbal frame structure has “negligible” inertia, 

6) angular rate of the gimbal frame is “negligible” com¬ 
pared to the rotor angular rate about its symmetry axis. 

Salient features of this VSCMG model are system adaptabil¬ 
ity (misalignment correction), and scalability. 

The ADCS mechatronics architecture comprises of a 
COTS smartphone running stock Android OS with self- 
contained ADCS routine, an ADK microcontroller board, 
ESC units along with the motor driven VSCMG arrays. The 
standalone mechatronics architecture performs the task of 
state sensing through embedded MEMS sensors, filtering, 
state estimation and implementation of feedback control law, 
to achieve the desired control objectives while maintaining 
active uplink/downlink with a remote ground control station. 

11. Attitude Dynamics of a Spacecraft with 
VSCMG ACTUATORS 

A. Spacecraft Attitude Kinematics 

Let R{t) denote the rotation matrix from the base body- 
fixed coordinate frame to an inertial coordinate frame. If n(f) 
is the total angular velocity of the base body with respect to 
the inertial frame and expressed in the base body frame, then 
the attitude kinematics of the base body is given by 

R{t) = R{t)n{t)^. (1) 

The spacecraft with a VSCMG has 5 rotational degrees of 
freedom, which are described by the variables; gimbal angle 
(a), rotor angle (9) and attitude R. 

B. Equations of motion 

The accurate dynamics model of a spacecraft with 
VSCMG derived in [7], [8], appropriate for model-based 
controller design, is briefly presented in this section. The 
gimbal axis (g) of the VSCMG is orthogonal to the rotor’s 
axis of rotation which passes through its CoM, if 

the rotor offset, cr = 0. A schematic CAD rendering of a 
VSCMG assembly, is shown in Eig. The VSCMG models 
accounts for rotor misalignments (a f 0) if any, which 
makes the system more adaptive and reliable. Total rotational 
kinetic energy of a spacecraft with a VSCMG is 
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Eig. 1; Exploded view of the VSCMG unit 
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Define the following angular momentum quantities, which 
depend on the kinetic energy Q; 

dT 

Total ^ n = — = (Jt + lT{t))^ + B{t)f, (7) 
dT T 

VSCMG ^p=— = B{t)'-VL +Jgr{t)f. (8) 

The equations of motion obtained from applying the 
Lagrange-d’Alembert principle can be written as follows; 

( 9 ) 
( 10 ) 

where, M.g{R{t)) is the gravity gradient moment on the 
spacecraft, expressed in the base body frame and r is the 
VSCMG control torque. The total angular momentum of 
the spacecraft can be expressed as the sum of the basebody 
momenta and the VSCMG momenta. 

Spacecraft basebody momenta 

n(f)= Ijt + iT{t))n{t) + Bit)j{t), (11) 

VSCMG momenta 


— =nxn + Mg{R{t)), 
dp dT 
di 


where, Jt = Jb - {rUg + mr)(pg)'^, 

Irit) = Jc{t) + Ic{t) - TTlrCrg^t)^Pg , (12) 

Jc{t) = Rg{t)JgRg{t)'^ + Rr (t) Jr Rr {t)'^ , (13) 
Ic{t) = -mr<j{pgp{t)^ -f CT(r7(f)^)^). (14) 


















Here, is the inertia of the spacecraft basebody; Jg and 
are the corresponding gimbal and rotor inertias, respectively. 
Total constant part of inertia in the base body frame is 
denoted by Jt whereas Ixit) is the total time-varying part 
of inertia in the base body frame. Jc{t) and Ic{t) are the 
offset-independent part of inertia and offset-dependent part 
of the VSCMG inertia in the spacecraft base body frame, re¬ 
spectively. Also, Rg and Rr denote the rotation matrix from 
the VSCMG gimbal-fixed and rotor-fixed coordinate frames 
to the spacecraft base body coordinate frame, respectively. 

The dynamics of the base body is therefore affected by 
the rotation rates a{t) and 9{t) of the VSCMG gimbal and 
rotor respectively, but is independent of the rotation angle 
0{t) of the rotor which is a cyclic variable for this system. 
Therefore, the conjugate momentum corresponding to 9{t), 
i.e. pe{t), is conserved if there is no torque input acting on 
the rotor axis. Note that the angular speed of the rotor 6{t) 
is not constant in this case due to the dynamical coupling 
between the different degrees of freedom of this spacecraft 
system.Therefore, to keep a constant rotor angular speed in 
the base spacecraft body frame as in constant speed SGCMG, 
a control torque has to be applied to the rotor. 


III. Spacecraft ADCS 

Spacecraft ADCS enables the spacecraft to realize its 
current states on SO(3)and can reorient the spacecraft to any 
desired final states. The raw IMU measurements from the 
smartphone are fused/filtered through the estimation scheme 
given in Section III-A Section III-B provides the control 
strategy for generic spacecraft attitude maneuvers. 


A. Estimator 


corresponding matrix of all these vectors expressed in the 
inertial frame as E = [ei 62 ei x 62 ]. Note that the matrix 
of the actual body vectors Ui corresponding to the inertial 
vectors e^, is given by ?7 = r"^E = [ui U2 ui x U2]. 

Proposition 1: A discrete-time filter that gives a second 
order numerical integrator for the filter in [15] is given as 
follows; 

4+1 =4exp(/i(H™i 

mwi+i = exp(-|Hf+i){("^-^3x3 - ^D)uJi+i 

h ^ - t (15) 

where, Wi+i =(m/ 3 x 3 + ^D) exp ( - ^nf)muji 

+ ^<i>'4"(4,C/r))^L.(4)| (16) 

is a discrete-time approximation to the angular velocity at 
time := ti + ^, where h is the time step, m is a 

positive scalar and D is a positive definite matrix, (4) = 
vex(Lj4 ~ RJLi) G and vex(-) : so(3) —>■ is the 

inverse of the map, Li = EiWi(UJ^)'^ G and 

(Rqj^o) G S0(3) X are initial estimated states. Further, 
W and <i)(ai) are chosen to satisfy the conditions in [15], and 

u\R,,ur) = R.ur,[E, - mr)w,), hd 

denotes the artificial potential term (Wahba’s cost function) 
which is generalized by applying the function T* on it. 


In [15], an estimation of rigid body attitude and angular 
velocity without any knowledge of the attitude dynamics 
model, is presented using the Lagrange-d’Alembert prin¬ 
ciple from variational mechanics. This variational observer 
requires at least two body-fixed sensors to measure inertially 
known and constant direction vectors as well as sensors to 
read the angular velocity. A first order discretized estima¬ 
tion scheme for computer implementation is also presented 
using discrete variational mechanics. Here, a second order 
symmetric Lie group variational integrator is introduced. 

In order to determine three-dimensional rigid body attitude 
instantaneously, three known inertial vectors are needed. This 
could be satisfied with just two vector measurements. In 
this case, the cross product of the two measured vectors is 
considered as a third measurement for applying the attitude 
estimation scheme. Let these vectors be denoted as u™ 
and u™, in the body-fixed frame. Denote the corresponding 
known inertial vectors as seen from the rigid body as ei and 
62 , and let the true vectors in the body frame be denoted 
Ui = R Ci for z = 1,2, where R is the rotation matrix from 
the body frame to the inertial frame. This rotation matrix 
provides a coordinate-free, global and unique description of 
the attitude of the rigid body. Define the matrix composed 
of all three measured vectors expressed in the body-fixed 
frame as column vectors, = [it™ it™ it™ x u™] and the 


B. Attitude Control using VSCMG array 
From o one can express, 

n = Hf, -f It, where H;, = A(f)f2 and u = BT. (18) 

The idea behind controlling the spacecraft attitude dynam¬ 
ics with the VSCMG is to ensure that the angular velocity 
of the spacecraft, which depends only on the term H^ in 
equation ( [TS] ), is controlled using the internal momentum it, 
which depends on the VSCMG rotation rates F. For simplic¬ 
ity, we consider the case where there is no gravity gradient 
moment, i.e., Mg{R{t)) = 0. Therefore, re-expressing the 
base-body attitude dynamics as: 

Hf, = H;, X H -f (it X H — li) = Hft X H -f Tcp, (19) 

where Tcp = it x H — it is the control torque generated by the 
“internal” momentum it from the VSCMGs. Note that from 
d and the attitude kinematics Q, the total momentum in 
an inertial frame, H/ = i?n = i?(nf,-|-ii) is conserved. There 
are several control problems of interest for a spacecraft with 
n VSCMGs. For most applications, the common spacecraft 
attitude maneuvers are slewing to rest, pointing , and attitude 
tracking maneuver. The results and discussion for a slew 
to rest attitude control scheme (de-tumbling maneuver) and 
a pointing maneuver for a spacecraft with n VSCMGs are 
presented in references [7], [ 8 ] and [16] respectively. 




IV. Mechatronics architecture 
A. Hardware 

For demonstrating the technological feasibility of the 
smartphone based ADCS using VSCMG, a lU cubesat with 
miniature VSCMG actuator array as shown in Fig. is 
considered, nevertheless the VSCMG system is scalable to 
any operational range from cubesat to space station. Each 
hermetically sealed VSCMG unit houses two brushless DC 
motor individually sized, controlled to drive the gimbal and 
rotor. An exploded of VSCMG showing its internal compo¬ 
nents is depicted in Fig. [T] The VSCMG array is formed by 
arranging 4 individual units in pyramid conhguration, where 
its skew angle is given by /3. The main advantage of using 
VSCMG as attitude actuator is to exploit its output torque 
to the input power ratio properties. 



Fig. 2: Spacecraft bus with smartphone and ADCS 

Modern smartphones are integrated with the following 
hardware sensors: 

1) Accelerometer 

2) Gyroscope 

3) Magnetometer 

4) GPS, 5) Barometer, 6) Light and proximity sensor etc. 
Out of these hardware sensors, (1), (2) and (4) are fused 
and used as an Inertial Measurement Unit (IMU) for attitude 
determination. Sensors (3) and (5) can be used for altitude 
measurements in UAVs. The estimation and controller rou¬ 
tine embedded in the on-board smartphone computes the 
required output signal and commands the VSCMG array 
to generate the corresponding actuator torque via Electronic 
Speed Control (ESC). The control signal from the on-board 
spamrtphone is sent to an microcontroller enabled with Ac¬ 
cessory Development Kit (ADK) compatibility (eg. Google 
ADK, Arduino MEGA ADK, lOIO, Microbridge etc.), where 
corresponding Pulse Position Modulation (PPM) signals are 
generated and sent to ESC to control each VSCMG param¬ 
eters, i.e. gimbal angle (a), gimbal rate (d(f)), rotor angle 
(9) and rotor rate (9{t)). The ESC units along with an Uni¬ 
versal Battery Elimination Circuit (UBEC) are independently 


calibrated and programmed to drive each gimbal and rotor 
motors of the VSCMG array. The microcontroller interfaces 
the smarphone with other peripherals including ADCS array 
and payloads; in addition to the system watchdog service. 
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Eig. 3: Hardware architecture of smartphone based ADCS 


B. Software 

1) Smartphone: Android is a Linux-based, open source 
mobile operating system that is very popular among software 
developers because of its powerful capabilities and open ar¬ 
chitecture. ADCS processes like sensor data acquisition, pre- 
hltering, sensor data fusion, estimation and control schemes 
are executed on-board the smarthphone’s processor. These 
functions are run as embedded native machine code compiled 
from C and/or C-H- source hies in the Android Application 
Packages (*.apk) using the Android Native Development Kit 
(NDK) toolset [17]. Android NDK-based spacecraft ADCS 
programming is especially useful because it enables carry¬ 
ing out CPU-intensive workloads in real-time. A low level 
architecture shown in Fig. is exploited for implementing 
the data acquisition from the accelerometer, gyroscope and 
magnetometer sensors. This is then used by the estimation 
scheme for attitude and angular velocity estimation through 
the Android Java application that interacts with native code 
using the Java Native Interface (JNI). External sensor break¬ 
outs and payloads are coupled to the application framework 
through the stable Linux subsystem I2C by writing the 
drivers in kernel and linking to the Hardware Abstraction 
Layer (HAL). 

2) Microcontroller: The Accessory Development Kit 
(ADK) is a microcontroller development board that adheres 
to the simple Open Accessory Standard Protocol created by 
Google as a reference implementation [18]. The ADK 2012 
is based on the Arduino open source electronics prototyping 
platform, with some hardware and software extensions that 
allow it to communicate with Android devices. The ADK 





































ANDROIDAPPLICATION 


1 


APPLICATION 

FRAMEWORK 

1 










HARDWARE ABSTRACTION LAYER 



KERN EL (DRIVERS) 




i 

k i 

k 




Accelerometer 

Gyroscope 

Magnetometer 

Camera 

USB 



Fig. 4; low level architecture 


microcontroller has embedded within it the system program 
written in the C++ IDE, in parallel to the Android side 
programming. In addition, a sophisticated power manage¬ 
ment scheme can be implemented in the microcontroller, 
as all the system components including the smartphone and 
VSCMG actuators are powered by external battery packs via 
the UBEC and microcontroller. 

V. HIE SIMULATION OF SPACECRAFT ATTITUDE 
CONTROL 

A. Testbed Hardware 

The spacecraft attitude simulator testbed uses hemispher¬ 
ical air-bearing to perform unrestricted yaw motion and re¬ 
stricted pitch and roll motion. The spacecraft bus is mounted 
on the air-bearing assembly and exhibits the same rotational 
degree of freedom (DoE) as that of the air-bearing. The 
spacecraft bus is actuated by a cluster of three VSCMG 
arranged in tetrahedron configuration. Each VSCMG is inde¬ 
pendently driven by a high torque stepper motor to drive the 
gimbal and a variable speed servo motor to drive the rotor. 
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Eig. 5: HIE testbed 


B. Experiments 

The estimation scheme given in Section III-A is imple¬ 
mented off-board on a remote PC using the sensor measure¬ 
ment acquired and transmitted by the on-board smartphone. 
The coordinates used as the inertial frame is ENU, which is a 


right-handed Cartesian frame formed by local east, north and 
up. The coordinates fixed to the COM of the cellphone with 
right direction of the screen as x, up direction as y and the 
direction out of screen as z is considered to be the body fixed 
frame. As mentioned in |III-A| at least two inertially known 
and constant directions are required in order to estimate 
the rigid body attitude. Given the sensors installed on the 
smartphone, one could use the accelerometer to measure 
the gravity direction and the magnetometer to find the 
geomagnetic field direction. The cross product of these two 
vectors is considered as the third vector. In order to find these 
directions, one could normalize the data from accelerometer 
and magnetometer in the case that the cellphone is aligned 
with the true geographical directions and the body fixed 
frame coincides with the ENU frame. Note that the direction 
read by the accelerometer shows the up direction, since an 
upward acceleration equal to g is applied to the phone in 
order to cancel the Earth’s gravity and keep the phone still. 
Therefore, the matrix of three inertially constant directions 
is expressed in ENU frame as 


E = [ei 62 61 X 62] = 


0 

0 

1 


0.0772 

0.6117 

-0.7873 


-0.9921 

0.1251 

0 


The three axis gyroscope also gives the angular velocity 
measurements. These three sensors produce measurement 
data at different frequencies. The filter’s time step is selected 
according to the fastest sensor, which is accelerometer here. 
At those time instants where some of the sensors readings 
are not available because of the difference in sampling 
frequencies, the last read value from that sensor is used. 

$(•) could be any function with the properties de¬ 
scribed in Section 2 of [15], but is selected to be $(a;) = x 
here. Eurther, W is selected based on the value of E, such 
that it satisfies the conditions in [15] as below; 


W = 


3.19 

1.51 

0 


1.51 0 
3.19 0 
0 2 


The inertia scalar gain is m = 0.5 and the dissipation matrix 
is selected as the following positive definite matrix: 

D = diag([12 13 14]'^). 


Sensors outputs usually contain considerable levels of 
noise that may harm the behavior of the nonlinear filter. A 
Butterworth pre-filter is implemented in order to reduce these 
high-frequency noises. Note that the true quantities would 
not contain high-frequency signals, since they are related to 
a rigid body motion. A symmetric discrete-time filter for the 
first order Butterworth pre-filter is implemented for filtering 
the measurement data as follows: 


(2 + h)xk+i = (2- h)xk + h{x^ + (20) 

where h is the time stepsize, x and x"* are the filtered and 
measured quantities, respectively, and the subscript k denotes 
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VI. Conclusion 




Fig. 7; Angular Velocity Estimation Error 


the fcth time stamp. The initial estimated states have the 
following initial estimation errors: 

Qo = expmso(3) (^(2-2 x [0.63 0.62 - 0.48]"^^ 

and Wo = [0.001 0.002 - 0.003]"^ rad/s. (21) 

In order to integrate the implicit set of equations in ( [TS] ) 
numerically, the first equation is solved at each sampling 
step, then the result for Ri+i is substituted in the second one. 
Using the Newton-Raphson method, the resulting equation 
is solved with respect to uJi+i iteratively. The root of this 
nonlinear equation with a specific accuracy along with the 
Ri+i is used for the next sampling time instant. This process 
is repeated to the end of the simulation time. Using the 
aforementioned quantities and the integration method, the 
simulation is carried out. 

C. Results 

Experimental results for the attitude estimation scheme, 
obtained from the experimental setup described in the previ¬ 
ous subsection, are presented here. These experiments were 
carried out on the HIE simulator testbed in the Spacecraft 
Guidance, Navigation and Control laboratory at NMSU’s 
MAE department. Experimental results from the control are 
not available yet since the VSCMG actuators have not been 
finalized and mounted on this testbed. The principal angle 
coiTesponding to the rigid body’s attitude estimation eiTor is 
depicted in Eigure [^ Estimation errors in the rigid body’s 
angular velocity components are shown in Eigure All the 
estimation errors are seen to converge to a neighborhood of 
(Q, w) = (/, 0), where the size of this neighborhood depends 
on the characteristics of the measurement noise. 


This paper presents a novel overall (hardware and soft¬ 
ware) architecture of a spacecraft attitude determination 
and control subsystem (ADCS), using a smartphone as the 
onboard computer. This architecture is being implemented 
using a HIE ground simulator for three-axis attitude motion 
simulation in the Spacecraft Guidance, Navigation and Con¬ 
trol laboratory at NMSU. Theoretical and numerical results 
for the attitude control and attitude estimation schemes 
that are part of this architecture, have appeared in recent 
publications. Both the attitude control scheme and the at¬ 
titude estimation scheme provide almost global asymptotic 
stability, and are robust to measurement noise and bounded 
disturbance inputs acting on the spacecraft. Experimental 
verification of the attitude estimation algorithm is presented 
here, and the experimental results show excellent agreement 
with the theoretical and numerical results on this algorithm 
that have appeared in recent publications. 
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